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Abstract. We consider a semiclassical 2x2 matrix Schrodinger oper- 
ator of the form P — — /i 2 Al2 + diag(Vi(a;), V2(x)) + hR(x, hD x ), where 
Vi, V2 are real-analytic, V2 admits a non degenerate minimum at 0, Vi is 
non trapping at energy 14(0) = 0, and R(x, hD x ) = (rj,k(x, hD x ))i<j ik < 2 
is a symmetric off-diagonal 2x2 matrix of first-order pseudodifferen- 
tial operators with analytic symbols. We also assume that Vi(0) > 0. 
Then, denoting by ei the first eigenvalue of —A + (V 2 "(0)a;, x) /2, and 
under some ellipticity condition on n,2 and additional generic geomet- 
ric assumptions, we show that the unique resonance pi of P such that 
pi = V 2 (0) + (ei + r 2 , 2 (0, 0))h + 0(h 2 ) (as h -> 0+) satisfies, 

Im Pl = ^h na+ ^-^ /2 f(hM^)e- 2S/h , 

where f(h,\n\) ~ Z)(,<m<£ fe,mh e (ln \) m is a symbol with / ,o > 0, 
S > is the so-called Agmon distance associated with the degenerate 
metric max(0, min(Vi, V2))dx 2 , between and {V\ < 0}, and n > 1, 
nr > are integers that depend on the geometry. 
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1. Introduction 

The theory of predissociation goes back to the very first years of quantum 
mechanics (see, e.g., \Kr\ \L& \ IZe ] I5t]). Rouhgly speaking, it describes the 
possibility for a molecule to dissociate spontaneously (after a sufficiently 
large time) into several sub-molecules, for energies below the crossing of the 
corresponding energy surfaces of the initial molecule and the final dissoci- 
ated state. From a physical point of view, one naturally expects that this 
(typically quantic) phenomenon occurs with extremely small (but non zero) 
probability. 

Despite the fact that statements concerning this problem are present in the 
physics literature for more than 70 years, the first mathematically rigorous 
result is due to M. Klein [Kl] in 1987, where an upper bound on the time 
of predissociation is given in the framework of the Born-Oppenheimer ap- 
proximation. More precisely, denoting by h the square root of the ratio of 
electronic to nuclear mass, M. Klein proves the existence of resonances p with 
real part below the crossing of the energy surfaces, and with exponentially 
small imaginary part, that is, 

|Im p\ = 0( e - 2 (l-^A) 

where S > is a geometric constant, e > is fixed arbitrarily, and the 
estimate holds uniformly as h goes to zero. 

In terms of probabilities, this result corresponds to give an upper bound on 
the transition probability between the initial molecule and the dissociated 
state. The purpose of this article is to obtain a more complete information 
on this quantity, and in particular a lower bound on it. More precisely, 
under suitable conditions, we prove that the imaginary part of the lowest 
resonance admits a complete asymptotic expansion of the type, 

lvi Pl = -h no+ ^ n ^/ 2 e- 2S l h V /^//(ln 1 )™, 

^— ' ' h 

0<m<£ 

in the sense that, for any N > 1, one has 

|Im pi + /l "0 + (l-nr)/2 e -2S/fc £ fi^h'Qn ^) m \ 

0<m<e<N 
= 0(/i n o + ( 1 -"r)/2+Af e -25/h^ 

where S > 0, uq > 1 and np > are all geometric constants, and where the 
leading coefficient /o,o is positive. 
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As it is well-known, the quantity Im p is closely related to the oscillatory 
behavior of the corresponding resonant state in the unbounded classically 
allowed region. Hence, the main issue will be to know sufficiently well this 
behavior. 

The strategy of the proof consists in starting from the WKB construction at 
the bottom of the well, and then trying to extend them as much as possible, 
and at least up to the classically allowed unbounded region. This is mainly 
the same strategy used in |HeSj2| for the study of shape resonances. 

However, from a technical point of view, several new problems are encoun- 
tered, because of the crossing of the electronic levels. 

The first one is that, at the crossing, the only reference on WKB con- 
structions is that of |Pe| . that has been done for a special type of matrix 
Schrddinger operators. In particular, it strongly uses the fact that only dif- 
ferential operators are involved. In our case, since our operator comes from 
a Born-Oppenheimer reduction, it is necessarily of pseudodifferential kind 
(see, e.g., jKMSW, MaSo]). As a consequence, our first step will consist in 
extending the method of [Pej to pseudodifferential operators. Unfortunately, 
this extension is far from being straightforward, and needs a specific formal 
calculus adapted to expressions involving the Weber functions. 

The second one is that, after having overcome the crossing, the symbols of 
the resulting WKB expansions do not anymore satisfy analytic estimates 
(usually needed in order to re-sum them, up to exponentially small error 
terms). In particular, this prevents us from using directly the constructions 
of |HeSj2| near the classically allowed unbounded region. Instead, we have 
to adapt the method of |FLM| that, without analyticity, allows us to extend 
the WKB constructions into the classically allowed unbounded region up to 
a distance of order (/iln |/i|) 2 / 3 from the barrier. This is not much, but this 
is enough for having a sufficient control, in this region, on the difference 
between the true solution and the WKB one. This is actually done by 
adapting the specific arguments of propagation introduced in l'l.Mj . where 
the propagation takes place in /i-dependent domains. 

In the next section, we describe in details the geometrical context and the 
assumptions. 
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In Section [3j we state our main result. 

Section [4] is devoted to the WKB constructions, starting from the well and 
proceeding away along some minimal geodesies, until crossing the boundary 
of the classically forbidden region. It is in this section that we develop 
a formal pseudodifferential calculus adapted to expressions involving the 
Weber functions. 

Next, in Section [5j we extend the well-known Agmon estimates to our pseu- 
dodifferential context. In this case, the main feature is that, since we can- 
not use general Lipschitz weight-functions, we replace them by /i-dependent 
smooth functions with bounded gradient, but with derivatives of higher or- 
der that can grow to infinity as h — > 0. 

In Section [6j we use these estimates in order to obtain a bound for the dif- 
ference between the WKB solutions and a solution of a modified problem, 
and this permits us to define an asymptotic solution in a whole neighbor- 
hood of the classically forbidden region (but only up to a distance of order 
(hln {hi) 1 / 3 from this region). 

Section[7]contains the a priori estimates and the propagation arguments that 
lead to a good control on the difference between the asymptotic solution and 
the actual one. 

Finally, Section [8] makes the link with the width of the resonance. Even if 
the idea is standard (practically an application of the Green formula: see, 
e.g., |HeSj2| ), here we have to be careful with the double problem that, on 
the one hand, we deal with pseudodifferential (not differential) operators, 
and, on the other hand, the magnitude of freedom outside the classically 
forbidden region is of order (hln \h\)^/ 3 as h — > 0. 



2. Geometrical Assumptions 
We consider the semiclassical 2x2 matrix Schrodinger operator, 

(2.1) P= ( ^ M + hR(x,hD x ) 

with, 

P . ■- - h 2 A + Vj(x) (j = 1,2), 
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where x = (xx, . . . , x n ) is the current variable in M. n (n > 1), h > denotes 
the semiclassical parameter, and R(x,hD x ) = (rj^(x, /i-Dx))i<j,fc<2 is a for- 
mally self-adjoint 2x2 matrix of first-order semiclassical pseudodifferential 
operators. 

Let us observe that this is typically the kind of operator one obtains in the 
Born-Oppenheimer approximation, after reduction to an effective Hamil- 
tonian (see [KMSW, IHIHo]). In that case, the quantity h 2 stands for the 
inverse of the mass of the nuclei. 

Assumption 1. The potentials V\ and V 2 are smooth and bounded on W 1 , 
and satisfy, 

(2.2) ^i(O) > and E = is a non-trapping energy for V\; 

(2.3) V\ has a strictly negative limit as \x\ —> oo; 

(2.4) V 2 > ; F 2 _1 (0) = {0} ; Hesst^O) > ; liminf V 2 > 0. 

\x\— >oo 




In particular we assume that V 2 has a unique non degenerate well at x = 0. 
This well is included in the island O which is the bounded open set defined 

as, 

(2.5) O = {x € M n ; V\{x) > 0}. 

Next, we define the sea as the set where V\{x) < 0. 

The fact that is a non-trapping energy for V\ means that, for any (x, £) £ 
p ] ^ 1 (0), one has | ex.ptH pi (x,!;)\ — > +oo as t — > oo, where Pi(x,^) := ^ 2 + 
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V\(x) is the symbol of Pi, and H pi := (V^pi, — V X P\) is the Hamilton field 
of p\. 



Such conditions (2.2)-(2.4) correspond to the situation of molecular predis- 
sociation as described in [Kl] . 

Since we plan to study the resonances of P near the energy level E = 0, we 
also assume, 

Assumption 2. The potentials V\ and V2 extend to bounded holomorphic 
functions near a complex sector of the form, Sr s := {x E C n ; |Re x\ > 
Rq , |Im x\ < <5|Re x\}, with Ro, 5 > 0. Moreover V\ tends to its limit at 00 
in this sector and Re V2 stays away from in this sector. 
Assumption3. The symbols rj^(x,^) for (j,k) = (1, 1), (1, 2), (2, 2) extend 
to holomorphic functions near, 

Sr ,6 ■= <Sflo,5 x U G C n ; |Im £| < max(5(Re x), ^Mo)}, 

Mo > sup mw.(Vi(x),V2(x)). 
ara<i ; /or any a £ N 2n , i/iey satisfy 

(2.6) d a r jjk (x,C) = 0((Re 0) uniformly on Sr 0jS . 



Now we define the cirque O as, 



(2.7) n = {x£ R n ; V 2 {x) < Vi(x)} 

Hence, the well is in the cirque and the cirque is in the island. 

We also consider the Agmon distance associated to the pseudo-metric 

(min(lq,F 2 )) + cte 2 . 

Such a metric is considered in Pettersson |Pe] . There are three places where 
this metric is not a standard one. 

First, near the well 0, but this case is well known. It has been treated 
by Helffer and Sjostrand [HeSjlj. It is also considered in Pettersson. The 
Agmon distance, 

(2.8) <p(x)=d(x,0) 

is smooth at 0. The point (x, £) = (0, 0) is a hyperbolic singular point of 
the Hamilton vector field H q2 , where 02 = £ 2 — ^2(^)5 an d the stable and 
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unstable manifold near this point are respectively the Lagrangian manifolds 

= V(/?(x)} and {£ = -Vip(x)}. 
Secondly on d£l, precisely at the points where V\ = Yi. This case has 
been also considered in Pettersson. At such a point, if one assume that 
VVi 7^ VV2, then, any geodesic which is tranversal to the hypersurface 
{Vi = V 2 ] is C 1 . 

Finally there is the boundary of the island dO, where V± = 0. This situation 
was considered in Helffer-Sjostrand |HeSj2| . We will follow them in a next 
assumption. 

Now we consider the distance from the well to the sea, that means to dO : 
(2.9) S = d(0,dO), 

Setting Bs = {x £ O, tp(x) < S} and denoting by Bs its closure, we also 
consider the set Bs H dO, that consists in the points of the boundary of the 
island that are joined to the well by a minimal d-geodesic included in the 
island. These points are called points of type 1 in |HeSj2| , and we denote 
by G the set of the minimal geodesies joining such a point to in O. 

We make the following assumption. 

Assumption 4. For all 7 G G, 7 intersects dVL at a finite number of 
points and the intersection is transversal at each of these points. Moreover, 
Wi ^ W 2 on 7 n dn. 

Let us recall that the assumption that is a non trapping energy for V\ 
implies that VVi 7^ on dO, and therefore that dO is a smooth hypersurface. 

We define the caustic set C as the closure of the set of points x £ O with 
if(x) = S + d(x,dO). In particular, the points of type 1 are in C. As in 
|HeSj2| we assume, 

Assumption 5. The points of type 1 form a submanifold V, and C has a 
contact of order exactly two with dO along T. 

We denote by nr the dimension of T. Moreover, for any 7 G G, we denote 
by := #(7 H d£l) the number of points where 7 crosses the boundary of 
the cirque, and we set, 

no := min iV 7 ; Gq := {7 G G ; -/V 7 = no}. 
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Then, we make an assumption that somehow insures that an interaction 
between the two Schrodinger operators does exist. 

Assumption 6. There exists at least one 7 6 Go for which one has the 
ellipticity condition 772(3;, iVy(x)) 7^ at every point x £ 7 n <9f2. 



3. Main Result 

Under the previous assumption we plan to study the resonances of the op- 



Op£(n,i) Op£(n, 2 ) 

Opf(rT^) Op£(r 2 , 2 ) 



erator P given in (2.1), where R(x,hD x ) is defined as 
R(x,hD x ) := 

where for any symbol a(x, £) we use the following quantizations, 

0p ^ (aMx)= (2^ / ^ x ~ y)i/ha MHv)dydt; 

Opf (a)n(x) = [ e«*-yW h a(y,Qu(y)dydt. 



(2irh) r 

In order to define the resonances we consider the distortion given as follows: 
Let F{x) e C°°(K n ,K n ) such that F{x) = for \x\ < R , F(x) = x for 
\x\ large enough. For 9 > small enough, we define the distorded operator 
Pg as the value at v = i9 of the extension to the complex of the operator 
U U PU~ 1 which is defined for v real small enough, and analytic in v, where 
we have set 

(3.1) U v cp{x) = det(l + vdF(x)) 1/2 <P{x + uF{x)). 

Since we have a pseudodifferential operator R(x, hD) the fact that U u PU~ l 
is analytic in v is not completely standard but can be done without problem 
(thanks to Assumption 3), and by using the Weyl Perturbation Theorem, 
one can also see that there exists £0 > such that for any 8 > small enough, 
the spectrum of Pq is discrete in [— Eo,£o] — i[O,eQ0]. The eigenvalues of Pq 
are called the resonances of P |Huj HeSj2 , [HeMaj . 



We will need another small parameter k > related to the semiclassical 
parameter h > 0, defined as, 

(3.2) k:=h1n-. 

In the sequel, we will study the resonances in the domain [sq, Ch] —i[0, Ck], 
where C > is arbitrarily large. In this case, we can adapt the WKB 
constructions near the well made in |HeSjl| , and show that these resonances 
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form a finite set {p±, ... , p m }, with asymptotic expansions as h — > 0, of the 
form, 

Pi ~ h ^2Pi4 hi/2 i 

where pjj S M and pjfi = ej + ^2 (0,0), ej being the j-th eigenvalue of the 
harmonic oscillator —A + (V^'(O)a;, x)/2 (actually, to be more precise, one 
must also assume that the arbitrarily large constant C does not coincide 
with one of the e^-'s). 

In this paper we are interested in the imaginary part of these resonances. 
We have, 

Theorem 3.1. Under Assumptions 1 to 6, the first resonance p\ of P is 
such that, 

Im pi = -h no+ ^~ n ^/ 2 f(h, In \)e~ 2S/k 
where f{h, In i) admits an asymptotic expansion of the form, 

/(Mni)~ Yl kmh'iln^r, (h^O), 

o<m<e 



with /oo > 0, and S > is defined in (2.9). 



Moreover the other resonances in [—£q, Ch] — i[0, Ck] verify 

Im Pj = 0{h^e~ 2S/h ), 
for some real j3j, uniformly as h — > 0. 



4. WKB CONSTRUCTIONS 

In this chapter, we fix some minimal d-geodesic 7 £ G, and we denote by 
. . . , x^ N ^ the sequence of points that constitute 7 n d£l, ordered from 
the closest to up to the closest to O (note that 7V 7 is necessarily an odd 
number). We also denote by 7^, 7^ 2 ', • • • , j( n ~i +1 ) the portions of 7\<9Q 
that are in-between and x^ 1 ), and x^ 2 \ x( N ~>' and O, respectively, 
in such a way that we have, 

7 = 7 W u U 7 (2) U • • • U {x^} U 7 ^ +1 \ 

where the union is disjoint (in particular, by convention we assume that 
£ 7 ( - 1 ^)- Moreover, we start by considering the first resonance p\ only. 
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4.1. In the cirque. As in |Pe| . the starting point of the construction con- 
sists in the WKB asymptotics given near the well x = by a method due 
to Helffer and Sjdstrand H eSjlj . More precisely, because of the matricial 
nature of the operator and the fact that p\ is elliptic above x = 0, one finds 
a formal solution w\ of Pw\ = p\W\ of the form, 

(4.D ».^")=(t 1 (L%) ) )^ w/ '- 



where (p is defined in (2.8), and cij (j = 1, 2) is a classical symbol of order 



in h, that is, a formal series in h of the form, 

(4.2) a,j(x, h) = ^ h k a jtk (x) 

k=0 

with a^fc smooth near (here no half-powers of h appear since we consider 
the first resonance p\ only). Moreover, ai is elliptic in the sense that a2,o 
never vanishes. Note that the generalization of the constructions of |HeSjl| 
to the case of pseudodifferential operators is done by the use of a so-called 
formal semiclassical pseudodifferential calculus, which in our case is based 
on the following result: 

Lemma 4.1. Let <p = <p(x) be a real bounded C°° function on IR n and let 
p = p(x,£) S S(l), that extends to a bounded function, holomorphic with 
respect to £ in a neighborhood of the set, 

{(x,f) G SuppV^ x <C n ; |Im£| < \V<p(x)\}. 

Then, denoting by Op^ the left (or standard) semiclassical quantization of 
symbols, the operator e^^ h Op^(p)e^^^ h is uniformly bounded on L 2 (M. n ), 
and for any a E Co°(IR n ) and N > 1, one has, 

(4.3) (e^"O v L Me-^ h a \(x;h) 



E - 



-V 1 d?p(x,iV(p(x))d« (a(y)e^ x ^ h ) +0{h N l 2 ), 
a\ \l J s ■ y \ j y=x 



\a\<N 

locally uniformly with respect to x, and uniformly with respect to h small 
enough. Here, &(x,y) := <p(x) — (p(y) — (x — y)Vip(x). 

The proof of this lemma is rather standard and we omit it (see e.g. |Malj ). 
Then, the construction can be performed by using the formal series given in 



(4.3) in order to define the formal action of R(x,hD x ) on w\. Afterwards, 



these constructions can be continued along the integral curves of the vector 
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field V^p2(x,i'Vip(x))D x = 2Vip(x).V x (that is, along the minimal geodesic 
of d starting at 0), as long as pi(x,iV<p(x)) does not vanish (that is, as 
long as these minimal geodesies stay inside the cirque f2). In that way, after 
resummation and multiplication by a cut-off function, we obtain a function 



w\ of the form (4.1), that satisfies, 

(4.4) P Wl - Pl W! = 0(h°°e-^ h ), 

locally uniformly in (J 7, where the union is taken over all the minimal d- 



geodesics 7 coming from the well and staying in Q.. In particular, (4.4) is 
satisfied in a neighborhood J\f± of 7W. 



4.2. At the boundary of the cirque. Now, we study the situation near 
the point x^> G dfl. By jPe], Theorem 2.14, we know that there exists a 
neighborhood V\ of x^ and two positive functions <£i,</?2 G C°°(Vi) such 
that, 

tpi = if on Vi n {Vi < V 2 }; 

ip 2 = tp on Vi n {V 2 < Vi}; 

|V^-(x)| 2 = ^(x) (j = l,2); 

tpi = (f2 and Vipi = Vif2 on Vi D 6T2; 

^2(x) — tpi(x) ~ <i(x,<9r2) 2 . 

Actually, ^2 is nothing but ^(0, x), where is the Agmon distance asso- 
ciated with the metric V2(x)dx 2 , and tpi is the phase function of the La- 
grangian manifold obtained as the flux-out of {(x, V^O^)) ; x S Vi n dQ} 
under the Hamilton flow of <Zi(x,£) := £ 2 — V\(x). 

Then, we set, 

(4.5) tp := 2 (Vi + V2) , 

and we consider the smooth function z{x) defined for x £ Vi by, 



(4.6) 



z(x) 2 = 2 (v?2(x) - ^i(x)) 
z(x) < on Vi n {V 2 < Vi}. 
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In order to extend the WKB construction (4.1 ) accross <9il near we still 



Vh 1 



follow |Pej and try a formal ansatz of the form, 

(4.7) 

w 2 {x; h) = Y,h k (a k {x, h)Y kfi + AS fe (:c, h)Y kj 



k>0 

where, 



(4-8) a k (x, h) = I ha ^) ) ; p k(x , h) = ( /W?' h l 
v ; fcV ' ; V «fe,2(x,/i) y ; V hp kj2 (x,h) 

a k> j and are formal symbols of the form, 
(4.9) ^ ^ ^(lnfcrV'^) 

;>0 m=0 

(with 7'' m smooth in ui\), and for any fc > and e G C, the function Yfc )£ is 
the so-called Weber function, defined by, 

(4-10) Y Ke {z) = d k e Y^ £ {z) 

where Yq >£ is the unique entire function with respect to e and z, solution of 
the Weber equation, 

(4.H) Y^ + ^-e-^)Y o , E = 

such that, for e > 0, one has, 

(4.12) Yo, £ (z) ~ ^e z2 ' A z e - 1 (z -> +oo). 



As it is shown in |Pej . Theorem 4.3, a resummation of (4.7) is possible up 
to an error of order 0(h°°e-^ h ) . 

Now, since cp is not C°° (but only C 1 ) near x^\ we need to find some 
generalization of lemma [4~Tj For technical reasons, in the rest of this section 
we prefer to work with the right semiclassical quantization of symbols, that 
we denote by Op^. 

For u > and g G C°°(R 2n ;M+), we denote by S uo (g(x,0) the set of 
(possibly /i-dependent) functions p G C°° (M 2n ) that extend to holomorphic 
functions with respect to £ in the strip, 

A„:=feOeK n xC n ; \1m£\<uo}, 

and such that, for all a G N 2n , one has, 

(4.13) d a p(x,0 = O(g(x,Re0), 
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uniformly with respect to (x, £) G A uo and h > small enough. We also 
denote by So(g) the analogous space of smooth symbols obtained by substi- 

We first show, 



tuting M 2n to and "smooth" to " holomorphic" 



Lemma 4.2. Let i> > 0, m G M, p = € 5 i/0 ((^) m ), and iet = 

be a reaJ bounded Lipschitz function on W 1 such that 

\\V<f)(x)\\ L ™ < v ] 

Let also a = a{x ; h) € C°°(R n ) be such that, for all a G N n , 

(hD x ) a a(x;h) = 0{e~^/ h ), 

uniformly with respect to h small enough and x G R n . Then, 

(Opf(p)a) (x;/ l ) = 0(e-^)/' 1 ) 

uniformly with respect to h small enough and x G M" - . 

Proof. - We write, 
(4.14) 

e <K*)/hOp%(p)a(x;h) = | e M£/M-^)/y y) ^)a{y ; 

and, following [Sj], we make the change of contour of integration in £, 

(4.15) R n 9^e + ^i^^, 

f - y\ 

where ||V0(a;)[|i,°° < v\ < Uq. We obtain, 

(4.16) e^ x)/h Op% (p)a{x ; h) 

with, 

6(x,y;h) = a(y ;h)e Wx) - ullx - yl)/h = o^)-<t>{y)-^~v\/hy 
Therefore, 

(4.17) e(x,y;h) = 0(e- s ^ h ), 
with 5 = v\ — ||V^>||l°° > 0. 



Then, in the case m < —n, the result follows immediately from (4.16 )-(4.17) 
(and standard estimates on oscillatory integrals). In the general case, we 
just write, 

(4.18) Op£(p) = Op^(p)(2^ - h 2 A x y k (2v - h 2 A x ) k , 



14 
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with k integer large enough (e.g. k = l + |[m]| + n), and, since Opt;(p)(2fo — 
h 2 A x )~ k is a semiclassical pseudo-differential operators with (/i-dependent) 
symbol in S Uo ((^) m ~ 2k ) C 5 , i/Q ((^) _n ~ 1 ), the result follows by applying the 
previous case with a replaced by (2uq — h 2 A x ) k a. □ 



Now, in view of defining a formal pseudodifferential calculus acting on ex- 



pressions such as (4.7), for j = 1, ...,n and x G u>\, we set, 



7 



Then, for any k > 0, we have (see |Pej formula (4.18)), 



(4-19) Mx):= [ dx J " ] GM 2 (R). 



(4.20) {hD Xj - iA^x)) ( \e-^' h 



v^fc / fc^fc-i i ( ^ 



A / 1 e ~^)/h_ 



If a and b are (scalar) formal symbols of the type (|4.9|), and fc G N, we set, 



(4.21) I k (a,b)(x ; /i) = a(x ; fc)y fc , ( + b(x ; ^ :( '' ) 



and we make act any diagonal matrix- valued function B{x) = diag(Bi(a;), Bz(x)) 
G .A/f 2 (IR) on I k (a, b)(x ; h) by setting, 

(4.22) B(x)I k (a, b)(x ; h) := h(B ia , B 2 b){x ; h). 

(this is possible since the Y k: o and F^i are linearly independent). Then, 
using ( |4.20 ), for all j = 1, ...,n we have, 

(4.23) (^. -iAj(x)) (l k (a,b)e-^ 

%(hD Xj a + VhbD X] z , hD Xj b) + 4_ x (0 , kVhaD Xj z)^j e~^ h . 

For A(x) = (Ai{x),...,A n {x)) G (.M 2 (IR)) n and a G N™, we also use the 
notation, 

(4.24) A(x) a = Ai(x) Ql ...A n (x)°" G 7W 2 (M), 

and we identify any £ G M n with (£1X2, ...,6»Zb) G (X 2 (M)) n . 
Then, we have, 

Lemma 4.3. Let Vq > sup xeuJi min (y/Vi(x), y/V 2 (x)) and m G K. Then, 
for any 5 = S(x,0 = diag(B 1 (x, 0, S 2 (x, £)) G A< 2 0M(£> m ))> fc > 0, a 
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and b in Cg°(u;i), and a G N n , one has, 
(4.25) 

Op£ (B(x, f ) (£ - iA(x)) a ) (l k (a, b) e -^ x ^ h ) = O (\\nh\ k h^ 2 e-^ x ^ h 



where the estimates holds uniformly for h small enough and x G W 1 . 

Proof. We prove it by induction on \a\, being careful to the fact that we 
do not have at disposal a symbolic calculus similar to that of the usual 
pseudodifferential operators (because of the special kind of action of the 
diagonal matrices on Ik(a,b), which actually does not commute with the 
oscillatory integrations). We first notice that, by |Pej - lemma 4.6, for (3 G N n 
and j G {0, 1}, one has, 

(4.26) (hD x f (Y kjj (^) e -^ h \ = O (\\nh\ k e-^' h ) . 



As a consequence, the result for a = follows directly from Lemma 4.2 
Now, assume it is true for \a\ < N (N G N fixed arbitrarily), and let 7 G N n , 
I7I = 1. Then, for \a\ < N, we have, 

(4.27) Op* {B(x, (£ - i A(x)) a+ ~<) 4(a, b)e~^l h 

= (2^ I ei{X ' V)mF ^ y ^ (Z- iA (.v)V ' h{a,b){jj)e-^/ h dyd^, 

with F a (y,£) = B(y,£) (£ — iA(y)) a . Now, assuming without loss of gener- 
ality that 7 = (1, 0, 0), and using the fact that, 

{ ie i(.*-vWh = -hD yi (e i( - x - y ^ h ), 

we obtain (denoting by F a> x and F a ^ the two diagonal coefficients of F a ), 
Op* {B(x, - iA(x)) a+ '<) I k (a, b)e-^/ h 

= (^rW / ei{X ~ y) * /h ( hD yi -iA x {y))I k {F a , x {y,Qa{y),F a , 2 {y^)b{y)) 



and therefore, by (4.23), 
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Op£ (B(x, (£ - iA(x)) a+ ~<) I k (a, b)e-^/ h 

= J^hyJ' el[X ~ y)i/h [h(hD yi F a , 1 a + VhF a , 2 bD yi z,hD yi F a>2 b) 

+4-1 (0, kVhF a)1 aD yi z)] e~'^ y)lh dydi 

= hOp%(F a )I k (D yi a,D yi b)e^/ h + hO^(D yi F a )I k (a,b)e^ h 
+VhOp%(amg(F a , 2 ,0))I k (bD yi z,0)e-^ h 
+VhOp% (diag(0, F a ,i))J fc _i(0, kaD yi z)e-^ h . 

Then, applying the induction hypothesis (and using the fact that D yi F a is 
a sum of terms of the type B'(y, ^)(r] — iA(y))" with |/3| = |a| — 1) this gives, 

Opf (B(x, (C - ^4(z)) a+7 ) I fc (a, b)e~*W h 

= O (\\xih\ k h l+ ^ + \hih\ k h l+ ^ + lln/i^/i 1 ^ + lln/i^" 1 ^ 1 ^) e~ v/h 



and the proof is complete. □ 



Now, for any smooth function / on M. n , we set, 

(4.28) f(A(x)) = diag f/(^M),/(^M) ) G ( - (R''..Vf 2 (C)). 



<9x 



Then, for any j> E Su ((O m ) an d for any TV > 1, Taylor's formula gives, 
p(x,0h= Yl ^p(x,iA(x))(a-iA(x)) a 

\a\<N 

+ Y, B a (x,0(ti-iA(x)) a , 

\a\=N+l 



(4.29) 



where I 2 is the 2x2 identity matrix, and the B a 7 s are in A4 2 (Sv ((£) m ))- In 



particular, using Lemma 4.3, for any a and b in Cq°(u>i), we obtain, 



Op* (p) (I fe (a,6)e-^) 
= E ^Opj?(^p(x,iA(x))(e-iA(x)) a )(/ fc (a,6)e 

H<iV 

+0(/i iV / 2 e-^ /l ) 
+0(/i 7V / 2 e-^ /l ), 



|a|<JV 
/3<a 
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and thus, writing down as before the corresponding oscillatory integral, in 
the same way we deduce, 

Op£(p) (I k (a,b)e-^ h ) 



P<a 



(4.30) +0{h N/2 e 



Now, for M E Z and f2 C 1" open, we consider the space of sequences of 
formal symbols, 

oo I 

S M (u 1 ):={a=(a k ) km ;a k (x,h) = ^ £ h l {hih)™^ (x) 5 

«=-Mm=0 

7i ra er( Wl )}. 

and, for a, 6 G S' a/ (cji), we set, 



(4.31) /(a, b)e~^ h := £ ^4(a fc , ^6 fe )e 



ip/h 

k>0 



Using (4.23), we see that, for j = 1, ...,n, the action of (hD Xj — iAj(x)) on 
such formal series satisfies, 

(4.32) (hD Xj - iAj(x)) I (a, b)e~^ h = /(L>, b))e~^ h , 

where Lj is the operator, 

(a,b) ^ 



L . , s M x S M -»■ S^ 1 x 5 M - X 



defined by, 



a{ = hD Xi a k + hb k D x ,z ; 



(4.33) b{ = hD Xj b k + (k + l)ha k+1 D Xj z, 

(k G N). In particular, using the notations L = (L±, L n ) and L a = 
L%\..L% n , for all a G N n , we have, 

(4.34) L a maps S M (^i) x S M (wi) into S M -M( Wl ) x S^H^wi). 

We also make naturally act any smooth diagonal .A/f 2(C)-valued function 
= diag(5i(x), J B 2 (x)) on S M x 5 M by setting, 

(4.35) B(a,b) = {B ia ,B 2 b), 
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and we define the formal action of a pseudodifferential operator with symbol 
p G Sv ((Q m )> on expressions of the type I(a, b)e~^^ h , by the formula, 

(4.36) Opf(p)(/(a,6)e-+"' N 
= E 

a6H™ 



I i{iA{x) + L) a -^A(xfd?p{x,iA(x))(a,b)) e~^ h . 



Then, in view of Lemma 4.3 and (4.30), we immediately obtain, 

Proposition 4.4. Let a,b G S M (u)\) and denote by I(a,b)e~^/ h any re- 
summation of I(a,b)e~^/ h up to a 0(h°°e~ v ^ h )-error term. Then, for any 
X G Co°( w i); the quantity Op^(p) ^x^( a > b)e~^l h ^j is a resummation of 
Op£ (p) (/(x^X^e - ^), up to a 0(/ i 00 e -^' l )-error term. 

In particular, the operator P naturally acts (up to 0(h°° e~ v ^ h )-error terms) 
on expressions of the type, 

(4.37, 

where ay = (a j:k ) k > and /3j = (Pj,k)k>o are in (j = 1,2). 

Writing down the equation Pw 2 = pxw 2 , setting, 

I 

i>0 m=0 

and the analog formula for /flj-jt, and identifying the coefficients of h l (lnh) m 

for < m < Z < 1, we find (denoting by v = ( ^ 1 ~?~ ^ iri < 1 ^ r i,2 

~ ~ ~ y /ir 2 ,i P2 + ^2,2 

the right-symbol of P), 



iV<^ 2 )a?;o + r i,2(x, iVtp 2 )a2§ + 



rV € pi(x,zV^i)(Vi 



(4.38) 

[^1(3;, iVy?i)Daj - i(V x .V^pi)(x,iV(fi) 
(4.39) +n,i(a;, *V^i) - pi] fi§ = 0; 



« = 0; 



r%>2(x,£V<p 2 )(V2) 



(4.40) 

(d^p 2 (x,iVip 2 )D x - i(V x .V^p 2 )(x,iV(p 2 



+ -((Ress^p2)(x,iVip 2 )Vz,V(ipi - ip 2 )) 



0.0 n 
a^ = 0; 



(4.41) 



+r 2 , 2 (x, iV</? 2 ) - Pi) "20 = °; 
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(Here we also have used the fact that p ~ J2k>i h k Pk &s h — > 0.) Identifying 
the other coefficients, one obtains a series of equations that (in a way similar 
to [Pej -section 4) can be solved in Vi (possibly after having shrunk it a little 
bit around x^), and in such a way that one also has, 

(4.42) w 2 -wi = 0{h co e- iplh ) locally uniformly in Vi n {V 2 < V{\. 



where w\ is defined in (4.1), and w±, w 2 are resummations of w\ and w 2 . 
Among other things, this implies, 



(4.43) 



a° o = a 2,o in Vi n {V 2 < V{\. 



Moreover, we see on M.39| and d4.4lb that (3^ (respectively a 2 ' ) are solu- 
tions of a differential equation of order 1 on each integral curve of the real 
vector field V(pi(y).V y (respectively Vip 2 (y).V y ). In particular, because of 
the ellipticity of 02,0, we deduce from ( 4.41| ) and (4.43) that we have, 



(4.44) 



a 2 'n never vanishes in Vi . 



Now, Assumption 6 implies that, if 7 G Go, then, 
(4.45) r h2 {x,iV<p 2 ) + on Vi- 



Since p\(y, iV(f 2 ) = Pi(y, iVtpi) = on uiDd^l, we deduce from (4.38) and 



(4.44) that, if 7 G Go, then /3^q does not vanish on oj\ n dft. As before, 
because of ( 4.39[ ) (and the fact that R(x, hD x ) is formally selfadjoint), this 
implies, 



(4.46) 



If 7 G G , then, ft 

'0 never vanishes in Vi . 



4.3. In the island, outside the cirque. Now, we look at what happens 
on 7^), and, at first, near x^ . Using the asymptotics of Y^^ e {z j \fh) given 
in |Pe]-section 4, one also finds that, in Vi n {V\ < V 2 }, w 2 can be formally 
identified with, 



(4.47) 



where b\, b 2 are symbols of the form, 



(4.48) 



b 3 {x-,h) =j2j2 hl ^ mbl r^ 

i>0 7Tt=0 



(j = 1,2), with b[: m G G°°(Vi n {Vi < V 2 }), in the sense that, for any 
resummations w 2 and W3 of w 2 and W3, one has, 

(4.49) w 2 - w 3 = O{h co e- Lpjh ) locally uniformly in Qn T+. 
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Moreover, one also has, 

(4.50) b{ =/3 1 ; 



which, by (4.46), shows that, when 7 6 Go, fri is elliptic in Vi n {V\ < V2}. 



Since P2(%, iV</?(x)) 7^ in {Vi < V2}, we can formally solve the equation 
P1V3 = P1W3, and we see again that b\ and 62 can be continued along the 
integral curves of Vip, as long as these curves stay inside {Vi < V2} and 
(pi does not develop caustics. In particular, they can be continued in a 
neighborhood M2 of 7^, and the continuation of b\ remains elliptic in ^2- 

Clearly, the previous steps can be repeated near x^ 2 \ 7^ 3 \ etc... (in the case 
i\L > 3), up to finally reach ■y( N ~< +1 \ obtaining in that way (after having 
pasted everything in a standard way by using a partition of unity) a function 
w(x, h), smooth on a neighborhood M("f) of 7 in O, and satisfying, 

(P-p 1 ) w = 0(h 0O e-^ h ), 

locally uniformly in M(j). Moreover, M (7) can be decomposed into, 

M(-y) = Mi U Vi U • • • U Vat 7 U M Nj +i, 

where, for all j, Vj is a neighborhood of and Mj is a neighborhood of 
7^, in such a way that, in each Mj, w admits a WKB asymptotics of the 
form, 

(4.5i) wfrh) -h^l h ZL/* )(x,h) ) ^ {x),h 



where af^ and are symbols of the same form as in fl4.48[), and a\ J > is 



elliptic if j is even, while 4 i} is elliptic if j is odd (in particular, a^ 7+1 ^ is 
elliptic). On the other hand, in each Vj, w can be representated by means 



of the Weber function, in a way similar to that of (4.7). 

4.4. At and after the boundary of the island. Let us denote by 

£ 7 6 7 n do, 

the point of type 1 where 7 touch the boundary of the island. When x £ 
7 H O is close enough to x 7 , we know from the previous subsection that the 
asymptotic solution w is of the form, 

(4.52) w(x; ft )~^(^"))e-^, 
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where 61,62 are smooth symbols on A/jv ' +1, of the same form as in (4.48), 
and 61 is elliptic. Moreover, as x approaches cc 7 , 61 and 62 (together with 
<p) develop singularities on some set C (called the caustic set). However, 
following an idea of |HeSj2| , we can represent ~e s l h w in the integral 
(Airy) form, 

(4.53) I[c ljC2 ](x,h) = h- 1 / 2 ^ ^ ( h^x\Z h h) ) e-^+^'^/^„, 

where we have used local Euclidean coordinates (x' , x n ) G M™^ 1 x M centered 
at 7n<90, such that V\(x) = — Cox ri + O(x 2 ) near this point. For x in O close 
to 7 n dO, the phase function £ n 1— >• x n £ n + g(x',£ n ) admits two real critical 
points that are close to 0. Then, choosing conveniently the x-dependent 
interval 7(3;), the steepest descent method at one of these points gives us 
the asymptotic expansion of I[c\,C2\. Comparing this with the symbols b\ 
and 62, one can determine c\ and C2 so that the asymptotic expansion of 
h 2 e b ' h "w coincides with that of I[c\, C2] in O. In particular, when 7 G Go, 
one finds that c\ remains elliptic near 0. 

At this point, since we did not assume any analyticity of the potentials near 
O, we have to follow the methods of jFLM] where a similar situation is 
considered. Indeed, following the constructions of |FLM] . Section 4 (that 
are made in the scalar case, but can be generalized without problem to our 
vectorial case), we see that there exists a constant 5 > such that, for any 
N > 1, one can construct a (vectorial) function wn, smooth on the set, 

(4.54) Wiv(7) := {\x - x y \ < e} n {dist(x, 6) < 2(Nk) 2/3 } 



with e > small enough (recall from (3.2) that k = \hlnh\), such that (see 
jFLMj . Proposit ions 4.5 and 4.6), 

• (P - pi)w N = 0{h 5N e- Kc * N l h ) uniformly in W N (-y); 

• For any a G ZIL, there exists m a > independent of N such that, 

uniformly in Wat (7); 



wn can be represented by an integral of the form (4.53) (with 7(x) = 
jn(x) depending on N) in all of Wat (7); 
w N = w in A/jv 7+ i n Wat (7); 

For any large enough L, there exist Cl > and 81 > 0, both 
independent of N such that, uniformly in Wjyi'y) H {dist(x,0) > 
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(Nk) 2/3 }, one has, 
(4.55) w N (x,h) 

^ ( L+[Nk/C L h] / ,c m 



\ 0<m<t 



£ h\\nh) m [ fl i N J:f ) + 0{h 5 ^ N + h L ) I e -Vvi-<-> '< 



as /i — >• 0, with fi™(x), /^(x) independent of /i, and of the form, 

(4.56) = (dist(x,C))- 3 ^ 2 - 1 / 4 /3^(x,dist(x,C)), 

(j = 1,2) where /J.'™ is smooth near (x 7 ,0), and /3i m (xj,0) 7^ in 
the case 7 G Go- 

Here, £>jv is a (complex-valued) C 1 function on Wat (7), smooth on Wn{i)\C, 
such that (see |FLM| . Lemma 4.1), 

• [p N = if + 0(h°°) uniformly in A/jv T +l H WaKt); 

• (V^jy) 2 = + C(/i°°) uniformly in Wtv(t); 

• There exists e(/t) = 0(h°°) real, such that, for x £ Wjy(7)\0, one 
has, 

(4.57) Re^jv(x) > S - e(h); 

• One has, 

Im V(/9at(x) = — i/jv(x)\/dist(x, C) Vdist(x,C) + C(dist(x, C)), 

uniformly with respect to h > small enough and x G VW(7)\0, 
with vn{x) > 5. 

The previous results show that we can extend w by taking wn in Wat (7), 
and we obtain in that way a function Wat smooth on A/"(7) U Wat (7), such 
that (P-pi)wjv = C(/i 5Af e- Re ^) uniformly in AA( 7 ) UWjv (7). Note that, 
thanks to Assumption 4, the number iV 7 is constant on each connected 
component of T. 

5. Agmon estimates 

5.1. Preliminaries. In order to perform Agmon estimates in the same 
spirit as in |HeSjl| , we need some preliminary results because of the fact that 
we have to deal with pseudodifferential operators (and not only Schrodinger 
operators). For this reason, we prefer to work with C°° weight functions (in- 
stead of Lipschitz ones) , and the idea is to take /i-dependent regularizations 
of Lipschitz weights. 
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At first, we need, 

Proposition 5.1. Let u > 0, m > 0, a = a(x,f) G S U0 ({£) 2m ). For h > 
smah enough, let also $>h £ C°°(IR n ) reaJ valued, such that, 

(5.1) sup|V$ h | < i/ , 
and, for any multi-index a S N n with \a\ > 2, 

(5.2) 5 a ^(x) = O (h 1 -^ 



uniformly for x G M n and h > small enough. Then, for any £ C R n 
with dist(S,M n \S) > 0, the operator e **/ fc Ar**/ A := e^OpJf (a) e -*"/ h 
satisfies, 

(5.3) Ue^Ae-^ullrn < Cill^Ac)™^^ 

uniformly for all h > small enough and u £ H rn (W a ). 

Proof. For it e Cg°(R n ), we write, 



and the property (5.1) shows that we can make the change of contour of 
integration given by, 

where ^(x, y) : = V<fr((l—t)x+ty)dt (in particular, one has: <fr(x) — &(y) = 
(x — y)^>(x,y)). Then, denoting by Op/j the semiclassical quantization of 
symbols depending on 3n variables (see e.g. |Ma2j Section 2.5), we obtain, 



and, using (5.2), we see that, for any a, /3, 7 6 Z+, we have 



(5.4) aztfqLfiM^ + Mfay))) =O{h-\^\0 m )- 

Then, the results is an easy consequence of the Calderon-Vaillancourt The- 
orem: see, e.g., [Ma2], exercise 2.10.15. □ 

We also need, 

Proposition 5.2. Let (ft and V be two bounded real-valued Lipschitz func- 
tions on W 1 , such that |V</>(x)| 2 < V(x) almost everywhere. Let also 
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Xi e C^(R n ; [0, 1]) supported in the ball {\x\ < 1}, such that J Xi(^)dx = 1. 
For any h > 0, we set Xh( x ) = h~ n x(x/h). Then, the smooth function, 

4>h ■= Xh * <P 

(where * stands for the standard convolution) satisfies, 

• 4>h = 4> + ®(h) uniformly for h > small enough and x € M. n ; 

• For all x G M n , one has \V(f) h (x)\ 2 < V{x) + h\\VV\\ L °o; 

• For all aeZf with \a\ > 1, one has d a (f> h = Oih 1 -^). 

The proof of this proposition is very standard and almost obvious, and we 
leave it to the reader. Observe that, in particular, <f>h satisfies the estimates 



(5.2). 



5.2. Agmon estimates. As a corollary of the two previous propositions, 
we have, 



Corollary 5.3. Let <j> and 4>h be as in Proposition 5.2, with V = min(Vi, V^)+ 
Then, for any u = (ui,u 2 ) G H 2 (R n ) H 2 (R n ), one has, 

2 

Re ^ h ' h Pu,e^ h u) > \\hV(e^/ h u)\\ 2 + J2((Vj ~ \V^ h \ 2 )e^/ h Uj , e^^uj) 

-C R /i(||e^ //l u|| 2 + ||/iV(e^ //l u)|| 2 ), 

where Cr > is a constant that depends on R(x,hD x ), xi and sup|V^| 
only. 

Proof. It is standard (and elementary) to show that, 

Re (e<^ //l (-/i 2 A + V^u^e^Uj) 

= UWCe^AijOU 2 + ((Vj - \Vc\> h \ 2 )e Mh Uj,e^l h Uj). 

Therefore, it is enough to estimate (e^ h ^ h R(x, hD x )u, e^ h / h u). Applying 



Proposition 5.1 , we see that the operator e^ h ^ h R(x,hD x )e ^ h / h {hD x ) 1 is 
uniformly bounded on L? . 



Moreover, since the constants appearing in the estimates (5.4) depend on 
a, a, and on the estimates on the (9^$'s only, we see that the norm of 
e ( t >h l h R( K x,hD x )e~^ h / h {hD x )~ l depends on r and on estimates on d@(xh * 
V0) = ip^Xh) * ( \/3\ < \a\) only. Since the latter depend on a, xi and 
sup|V</>| only, the result follows. □ 
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6. Global asymptotic solution 

The constructions of Section[4]can be done in a neighborhood of any minimal 
geodesic 7 G G, and give rise (after having pasted them together with a 
partition of unity) to an asymptotic solution (still denoted by w^v) on a 
neighborhood of lL e g 7- Now, we plan to extend this solution to a whole 
(/i-dependent) neighborhood of {V\ > 0}, by using a modified selfadjoint 
operator with discrete spectrum near 0. 

At first, we fix £0 > sufficiently small, and a cut-off function xo £ 
Cg°(0; [0,1]) such that, 

Xo(x) = 1 if Vi(x) > 2e ; xo(x) = if Vi(x) < e 0) 

and we set, 

(6.1) Vi :=x Vi + eo(l-Xo). 

In particular, Vi coincides with V\ on {Vi > 2e$}, and we have Vi > £o ev ~ 
erywhere. Then, we define Pi := — h A+Vi, and we consider the selfadjoint 
operator, 

(6.2) ?= ( P Q P 2 )+hR(x,hD x ). 

By construction, for all C > and h small enough, the spectrum of P is 
discrete in [— Ch,Ch], and a straightforward adaptation of the arguments 
used in |HeSjl| shows that its first eigenvalue E\ admits the same asymp- 



totics as p\ as h — > 0+. We denote by v its first normalized eigenfunction, 
and by Mq C {V\ > 2eq} some fix neighborhood of \J j& g n{V\ > 2eq} where 
the asymptotic solution wa? is well defined. We have, 

Proposition 6.1. There exists 9o£t independent of h, such that, for any 
compact subset K of Mo, and for any a S Z", one has, 

\\e^ h d a (e ido v - h^w N )\\ K = 0{h°°). 

Proof. The existence of #o such that d a (e ie °v -/i?Wjv) = 0(h°°) uniformly 
near 0, is a consequence of [HeSjlj, Proposition 2.5, and standard Sobolev 
estimates. Let \ £ C x> (Mq; [0, 1]), with x = 1 in a neighborhood of KL){0}. 



Following |HcSjf , IPej . we plan to apply Corollary 5.3 to u := x{ eie ° w 



/i4wjv), with a suitable weight function cf>. Let us first observe that, using 



26 ALAIN GRIGIS 1 & ANDRE MARTINEZ 2 

Corollary |5.3[ for any e > 0, one has, 

(6.3) \\e^- £ ^ h (hD x M m = 0(1). 

where <p(x) > <p(x) is the Agmon distance associated with min(Vi,V2) be- 
tween and x. Now, for C > 1 arbitrarily large, we define, 

:= min(0i,^ 2 ), 



where, 

4>l(x) :-- 



f{x) - Ch ln(yj(s) //i) if cp(x) > Ch; 
if{x)-ChhxC ii(p\x)<Ch, 



( inf {l-2e)(<p(y) + d{y,x)) if x £ suppx; 
^ 2 (x) := ^ X(2/)^1 

[ (1 — 2e)<^(:e) if x ^ suppx- 

Here, e > is taken sufficiently small in order to have <j>2(x) > (p(x) when 

x € K. Then, eft is Lipschitz continuous, and one has <f> = fa on X, and 

= 02 on M n \{x = 1}- Moreover, one sees as in |Pej (proof of Theorem 

5.5) that, if we set V := min(Vi, V2), 4> satisfies, 

\V<j)\ 2 = V in {ip < Ch}; 

|Vc/>| 2 < V - 5 Ch in {ip > Ch}, 

where 5q = inix&uppx; x^=o(V (x) / <p(x)) > 0. As a consequence, by Proposi- 



tion 5J2 the regularized fa of 4> satisfies, 

\Vfa,\ 2 < V + h\\V\\ L oo in {ip < Ch}; 
IV^I 2 <V— {5 C - \\V\\ L oo)h in {cp > Ch}. 
Then, choosing C sufficiently large, and setting u := x(e* e °v — h^Wx), we 



see that Corollary 5.3 implies, 

(6.4) \\hV{e^ h ' h u)f ^C'h\\e^ h uf { ^ Ch} < {e^\P - E ± )u, e^u), 

with C = C'(C) arbitrarily large. Moreover, if % € Cq°(A/o) is such that 
XX = X, we have, 

(P - E x )u = [P, X ]xu + 0{h°°e-^ h ), 

and since fa = (1 — 2e)</? + 0(h) on suppVx, min supp v x f =■ 5i > 0, and, 
by Proposition 5.1, the operator e^ h ^ h [R, x] e ~^ h ^ h is uniformly bounded, we 



obtain (using also (6.3)), 

( e ^A(P _ Ex)u, e^l h u) = O (\\e {1 ~ e) * /h (hD x )xu\\ 2 nppVx + h\\e^/ h u\\ 2 

= o(e~ 2e ^ h + h\\e^ h u\\ 2 ). 
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Inserting this estimate into (6.4), and taking C sufficiently large, this permits 
us to obtain, 

||fcV(e* h / h u)|| 2 + h\\e^ h/h u\\ 2 = 0(e- 2£5l/h + ||e* fc / h u||^< Ch} ). 

In particular, since 0^ = <f>\ + C(ft') on if, and 4>h = (1 — 2e)y3 < C7i on 
< Ch}, 

\\h c ^ c e^ h hVu\\ 2 K + [^V^N*^ = 0{e~ 2 ^l h + \\u\\\^ Ch} ). 
Therefore, 

He^Vull^ + \\e^ h u\\ 2 K = 0(h°°), 
and the result follows by standard Sobolev estimates. □ 



Now, following [FLMJ, Section 4.3, we observe that, if eo has been taken 
small enough, the asymptotic solution wjv is 0{h SN e~ s / h ) uniformly on the 
set, 

{dist(x, (J 7 ) > e } n {Vi < 2e } n |J AA( 7 ) U W N {l) 
7SG \ 7 eG 



Moreover, by (6.3), the same is true for v on {dist(x, 7) > eo} n {V\ < 

7£G 



2eq}. Therefore, using also Proposition 6.1, we can paste together e !e °v and 



h ™/ 4 W7v in order to obtain a function u^r that satisfies to the properties of 
the following proposition (see also |FLM| . Proposition 4.6): 

Proposition 6.2. There exists a function un, smooth on On '■= {dist(x, O)} < 
2(Nk) 2 / 3 , such that, 

5N -Re <p N /h\. 



(P-p)ujv = O(h oiy e 



d a u N = 0{h- ma e 



m a -Re tfN/h} 



uniformly on On, where ipN is as in (4.57). Moreover, un can be written 



as in (4.55) in [j Wjv( 7 ) n {dist(x,0) > (Nk) 2/3 } (with ,flf m (x 7 ,0) ^ 0), 



7SG 



while u N is 0{h 5N e~ Ke w/h) away from [j W N {l) D {x (£ 6}. 

■yeG 

7. Comparison between asymptotic and true solution 



7.1. A priori estimates. In the same spirit as in |FLM| . Theorem 2.2, we 
start with an a priori estimate for the resonant state of P. From now on, 
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we denote by u the out-going solution of 

(7.1) Pu = Pl u, 

normalized in the following way: we fix some analytic distorted space (also 
known, more recently, as a Perfectly Matched Layer), of the form, 

(7.2) R% := {x + i6F{x); x el"}, 

where F £ C°°(IR n ; W 1 ), F(x) = if \x\ < R , F(x) = x for \x\ large enough, 
and where 9 > is sufficiently small, and may also tend to with h, but not 
too rapidly (here, we take = h\mh\ = k). Then, by definition, the fact 



that p\ is a resonance of P means that Equation (7.1) admits a solution in 
L 2 (Mg), and here we take u in such a way that, 

(7.3) II U Hl 2 (R™) = 1- 

As before, d stands for the Agmon distance associated with the pseudometric 
min(Vi, V 2 )+dx 2 , and we denote by B d (S) := {x € R n ) ; d(0,x) < S} the 
corresponding open ball of radius S = d(0, dO). Then, we first have, 

Proposition 7.1. For any compact subset K C M. n , there exists Nk > 
such that, 

\\e s ^ h u\\ H1{K) = 0(h- N «), 
uniformly as h — > 0, where s(x) = ip(x) ifx G B^S) and s(x) = S otherwise. 

Proof. The proof if very similar to that of |FLM| Theorem 2.2, with the 
only difference that here we have to deal with pseudodifferentilal operators, 
forbidding us to use Dirichlet realizations and non smooth weight functions. 



Instead, we modify V± in a way similar to (6.1 ), and we regularize the weights 



as in Proposition 5.2 



We consider a (/i-dependent) cutoff function x such that, 

X (x) = 1 if Vi{x) > 2k 2 / 3 ; X (x) = if Vi(a?) < k 2 ^ ; d a X = 0(£T 2 l a l/ 3 ), 
and we set, 

Vi :=m + 3 {l-x) ; A :=-/i 2 A + t>i; 
(7.4) P= ( F Q p 2 )+hR(x,hD x ). 

We denote by E the first eigenvalue of P, and by v its first normalized eigen- 
function. Moreover, we consider the Agmon distance d associated with the 
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pseudometric ^min(Vi, — E\ dx 2 , and we set 0(x) := d(0,x). Then, 
the same proof as in [FLMJ, Lemma 3.1, shows the existence of a constant 
C\ > such that, 

(7.5) s(x) - dk < <p(x) < if(x) (x e M n ). 

Moreover, an adaptation of the proof of |FLM| . Lemma 3.2 (obtained by 



using Proposition 5.2 in order to regularize the Lipschitz weight) gives, 
(7.6) \\e^ h v\\ H Hun) = 0(h~ N °), 

for some Nq > 0. Then, the result follows by considering the function \v, 



and by observing that, thanks to (7.6), one has (see [FLMj, Lemma 3.3 and 
Formula (3.20)), 

||X* - ^ J(* ~ dz\\ Hl = 0(h- N ^e- s / h ). 

Here, 7 is the oriented complex circle {z G C ; \z — E\ = /i 2 }, and Pg is a 
convenient distortion of P. The previous estimate actually shows that the 
distorted ug of u coincides - up to 0(h~ Nl e~ s / h ) - with fxxv, where fi is a 
complex constant satisfying = 1 + 0(e~ s ^ h ), for some 5 > 0. □ 

Remark 7.2. The previous proof also gives a global estimate on ug, 

\\e s ^/ h ug\\ H1{Rn) = 0(h- N ^, 
for some constant N[ > 0. See |FLM] , Lemma 3.3 and Formula (3.20). 



Now, we plan to give an even better a priori estimate on the difference u— ujy 
near the boundary of the island. Here again, we follow the arguments given 
in [FLMj, Section 5. For any N > 1, we set, 

U N := {xeR"; dist(x, dO) < 2(Nk) 2/3 }. 

We have (see [FLMj , Propositions 5.1 and 5.2), 

Proposition 7.3. There exists N\ > and C > 1 such that, for any N > 1 
large enough, one has, 

\\u-u C N\\HHU N )<h- N2 e- S/h . 
Proof. We just recall the main lines of the proof in [FLMj . At first, thanks 



to Proposition 7.1 and the particular form of ucn, we immediately see that 
the estimate is true on the set {(f(x) > S — 2k}. Then, we take a cutoff 
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function x G C^{ip(x) < S — k) such that \ = 1 on {</?(x) > 5* — 2/c}, and 
d a x = 0(h~ Na ) for some A Q > 0. We also consider the Lipschitz weight, 

<p N (x) = min (ip{x) + dNk + k(S - v?(x)) 1/3 , S + (1 - fc 1/3 )(<S - <p(x)))\ , 



and, by using Propositions |7.1| and |6.2[ we see that if C is large enough, we 
have, 

\\ e ^/h {p _ pi) ~ (u _ UcN )\\ L2{Rn) = (h- M i), 
for some M\ > independent of N. Then, regularizing </>jv as in Proposition 



5.2, we can perform Agmon estimates as in the proof of [FLM], Proposition 

5.1, and we find, 

||e^x(u-ucAr)|U 2{R n ) =0(h-^), 
for some M\ > independent of N, and the result follows. □ 

7.2. Propagation. Now, we plan to prove (see [FLM], Proposition 6.1), 

Theorem 7.4. For any L > and for any a G Z" , there exists Nl j01 > 1 
such that, for any iV > Ni, a , one has, 

(7.7) ^(u-ucw)(i,fc) = 0(h L e - s/h ) as h^O, 

uniformly in Un- 

Proof. As in [FLM], the proof relies on three different types of microlocal 
propagation arguments. We fix some x G dO, and we define the Fourier- 
Bors-Iagolnitzer transform T (see, e.g., Sj IMa2j ) as, 



Tu(x,S;h):= [ e i( - x - y ^ h -^^ 2h u(y)dy. 
1) Standard C°° propagation 

Since u is outgoing (that is, it becomes L 2 when restricted to distorted 



space or Perfectly Matched Layer defined in (7.2)), one can see as in [FLMJ, 



Lemma 6.2 that, if to > is large enough, then, one has, 

Tu(x,0 = 0(/i 00 e- 5//l ), 



uniformly near exp(— toH pi )(x, 0). Moreover, by Proposition 7.1 we know 
that e s / h M remains 0(h~ N °) (for some Ao > 0) on a neighborhood of the 
x-projection of {exp(— tH pi )(x, 0) ; < t < to}. 
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Then, the standard C°° propagation of the Frequency Set for the solution 
to a real principal type operator (see, e.g., |Ma2j ) shows that the previous 
estimate remain valid near exp(— tH pi )(x, 0) for any t > 0. 

2) Non standard propagation in h-dependent domains 

Thanks to the previous result, we can concentrate our attention to a suf- 
ficiently small neighborhood of x. As before, we choose local Euclidean 
coordinates (x',x n ) G x M centered at x, such that V\{x) = —Cox n + 

0(\x — x\ 2 ). We also set u.n '■= (iV/c) -1 / 3 , and we considered the modified 
Fourier-Bors-Iagolnitzer transform Tjv defined by, 

(7.8) T N u(x,£;h):= [ e i ^- y ^ /h - ix '- y ' )2/2h - flNix "- y " )2/2h u{y)dy. 

Then, using the previous result it is elementary to show that (see [FLM], 
Lemma 6.3), for any (fixed) t > small enough, one has, 

T N l Kl VL{x,t) = 0{h co e- s l h ), 

uniformly near exp(— tH Pl )(x, 0),. Here K\ is of the form K\ = K\Bd(S), 
where K is any compact neighborhood of the closure of O. The interest 
of the latter property is that, as shown in [FLM], it can be propagated up 
/i-dependent times t of order (Nk) 1 ^ 3 . More precisely, setting, 

exptH pi (x, 0) = {x'{t),x n (t); ?(t),t n (t)) (t € R), 

we have (see [FLM] . Lemma 6.4), 

Lemma 7.5. There exists Sq > such that, for any 5 £ (0, So], for all N > 1 
large enough, and for tw,6 '■= 5 _1 (A^A;) 1 / 3 , one has, 

T N l Kl u = 0{h SN e- s/h ) uniformly in W(t N , h), 

where, 

W S (N, h) := {\x n - x n (-t N , 5 )\ < 5(Nk) 2 / 3 , |£„ - U~tN,s)\ < 5(Nk) 1 ^ 3 , 
W ~ x'Hn, S )\ < d(Nk) 1 / 3 , \£ - t'(-t N>s )\ < 5(Nk)^ 3 }. 

Proof. The proof is based on the refined exponential weighted estimates (in 
the same spirit as in [Ma2j) given in [FLMJ, Proposition 8.3 , that we apply 
here to the operator Pi. Since the proof is very similar to that of [FLMJ, 
Lemma 6.4, we omit the details. □ 
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On the other hand, using the explicit form of ucn given in (4.55), one also 
sees that, for any L large enough, there exists 5l > such that, for any 
N > 1, one has (see |FLMj . Lemma 6.7), 

T N l Kl u CN = 0{{h &LN + h L )e- s/h ) uniformly in W S (N, h). 

In particular, taking TV = L/Sl with L » 1, we obtain a sequence N = Ni 
along which, 

T N l Kl u CN = 0(h SLN e- s/h ) uniformly in W 5 (N,h), 
and with both N and 5lN arbitrarily large. 
As a consequence, along the same sequence, we also obtain, 

T N l Kl (u - u CN ) = 0(h 5 'L N e- s/h ) uniformly in W s (N,h), 
with 5' L = min(<5, 5l). 

Moreover, we see that, when y £ Un n Bd{S) and x G H x Ws(N, h) (where 
LTa; stands for the natural projection onto the x-space), we have, 

m( x n ~ Vnf + s(x) - S > C s Nk, 
with C$ > constant (and actually, C$ — > oo as 5 — > 0). Therefore, using 



Proposition 7.3 and the expression ( |7.8| ) of T^r, we also obtain, 

T N l UNnBd{s) (u - u CN ) = 0{h m e- s ' h ) uniformly in W S (N, h). 
as a consequence, if we set, 

(7 - 9) xn{x) = Xo » 

where the function xo £ Co°(^+; [0, 1]) verifies Xo = 1 i n a sufficiently large 
neighborhood of 0, and is fixed in such a way that xn(x) = 1 in {\x n — x n \ < 
\x n {-t N ) ~ x n \ + 25{Nk) 2 / 3 ; \x' - x'\ < \x'{-t N ) - x'\ + 25(Nk) 1 / 2 } (here, 
t jv and 5 are those of Lemma |7.5[ ), then, the function, 

vn ■= XNe S/h (u - u C n), 

is such that, 

(7.10) T N v N = 0(h s 'L N e~ s/h ) uniformly in W S (N, h). 

Moreover, we have (see jFLMj . Section 6.2), 

(P - Pl )v N = [P, X N]e S/h (u - u CN ) + 0(h 5N ), 



RESONANCE WIDTHS FOR THE MOLECULAR PREDISSOCIATION 33 

and thus, on {dw(x, suppVx/v) > e} x M n (where e > is fixed small 
enough and dw is the distance associated with the metric {Nk)~ 1 {dx 1 ) 2 + 

T N (P- Pl )v N = 0(h s ' N ), 

for some 5' = 5'(e) > 0. 

3) (Almost) standard analytic propagation 

Although we are in a region where no analytic assumption is made, a re- 
scaling of the problem makes appear estimates similar to those encountered 
in the analytic context. Indeed, setting, 

~ h =~ hN: =m = { m 4 

and performing the change of variable (still working in the same coordinates 
for which x = 0), 

x t-> x = (x, x n ) := {{Nk)- 1/2 x, {Nky 2/3 x n ); 



we see that the estimate (7.10) implies (see [FLMj . Formula (6.43)), 
uniformly in the tubular domain, 

W(h):={ IXn-Xni-S- 1 )] <S, ||n-Jn(-Ol<^ 

\x' - ^c-*- 1 )! < 5{Nk)-h , |f - e'C-^" 1 )! < s(Nk)-*y, 

where, 

v N (x) := (Nk^+^VNdNk^x' , (Nk)lx n ); 

(x(t),^(t)) :=exptfl~ (0,0); 

pi(x,0 := (iVfc^^p+^ + w^^X); 

WifoTi) := (7VA;)- 2 / 3 Fi((iVA ; ) 1 / 2 ^ 5 (iVA;) 2 / 3 J n ) - (Nk)~ 2 / 3 Pl . 
Moreover, setting, 

p ■= -(Nk)^ 3 h 2 A^ - h 2 dl n + Wi{x), 
then, for any N > 1 large enough, we also have, 

TPv N {x,lh N ) = 0{e- 5 ' ' 2 ~ h "), 
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uniformly with respectto h > small enough and (x, £) 6 M 2n verifying 
^(((iVfc) 1 ^, (jVA;) 2/3 5?„),suppVxi\r) > e. 



Finally, by Proposition 7.1 and Proposition 7.3 we have the a priori esti- 
mate, 



\VN 



for some iVj > independent of A, and we observe that, for iV = L/5l, one 
has Ni/(5 L N) -> as L -> +oo. 

At this point, a small refinement of the propagation of the microsupport (see 
jFLMj . Proposit ion 6.8) gives the existence of a constant 5i > independent 
of L such that, for all L large enough and N = L/Sl, one has, 

(7.12) Tv N (x,lh) = 0(e- 5 ^~ h ), 

uniformly in V{5 X ) = {x; \x\ < St} x (iVJfe)s|C'| + ||^| < <$i}. 

Then, using an ellipticity property of pi away from {^; (A£:)6|£'| + |£ n | < 
Si}, and reconstructing uat from Tvn, one finally finds, 

\\vN\\ H ™ { m< S2 )=0(e- s * s ^), 

with m > arbitrary, 52 > independent of L, N = L/Sl, and L arbitrarily 
large. Therefore, turning back to the original coordinates x and parameter 
h, and making x vary on all of dO, Theorem 7.4 follows. □ 



8. ASYMPTOTICS OF THE WIDTH 

As before, we denote by Pg the distorted operator obtained from P by means 



of a complex distortion as in (7.2), with Rq sufficiently large in order to have 
O C {\x\ < Rq/2}. We also denote by ug the corresponding distorted state 
obtain from u by applying the same distortion (see, e.g., [FLM] for more 
details). 

Let £ Cg°([0, 2); [0, 1]) with ^ = 1 near [0, 1], and set, 

/ dist(x,Q) ^ 

where, as before, N = L/Sl with L > 1 arbitrarily large. 
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Then, since ipN^ = 4>N~ue, Pe^e = Piug, and ipNPeipN~uo = ipNPipNU, we 
have, 

Im /jiII^atuH 2 = Im (ip N P e u e , ip N u) = Im ([ip N , P e ]u e ,ip N u), 

and thus, 
(8.1) 

(2/i 2 (VV'iv)Vu + /i 2 (A^)u,^u) + h([i; N ,Re]ug,^ N u) 
Im pi = Im — r-5 . 

Moreover, we know that ||t/>atu|| = l + 0(e~ s / h ) with 5 > and, by Theorem 
on supp ipNi we can replace u by ucn, up to an error 0{h L e~ s l h ). Also 



7.4 



using Proposition |7.1[ we deduce, 

(8.2) Im Pl = Im (2h 2 (V^ N )Vu CN + /i 2 (A^)u, Vtvuctv) 

+Mhfor, Re]u e , tPnu) + O(h L - N0 )e~ 2S / h , 
for some fix Nq > independent of L. 

Now, we introduce ipo G Cg°((l, 2); [0, 1]) with ijj = 1 near suppVf/'o, and 
we set i> N {x) = ^0 {j^0w) ■ 

Lemma 8.1. One has, 

(8.3) {[^ N ,Re\M e ^ N M) = & N [i> N ,R$ N u,$ N u) + 0(h°°e~ 2S / h ). 

Proof. Thanks to Assumption 3, we can make in [ipN,Rg] the (complex) 
change of contour of integration, 

x-y 



y/(x - y) 2 + h 2 

We obtain, 

ty N ,R e ]n e {x) = -^- J e^-y^ h -^l h )^ N {x)-^ N (yW e u e {y)dydi 
with, 

$ := yWo^t^L^ ; d^d^x, y, = (0). 
y / (x — y) z + h^ 

By construction, on the set, 

Aat := supp(Vw(x) - ^Ar(y)) n {y>jvO*0 / 1 or ^(y) / 1}, 

we have |x — y\ > c(Nk) 2 ^ 3 for some constant c > 0. As a consequence, on 
this set, the quantity \x — y\j \J{x — y) 2 + h 2 tends to 1 uniformly as h — > 0. 
Moreover, still on this set, we have either s(x) = S or s(y) = 5, and since 



86 
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\s(x) — s(y)\ < p\x — y\ with < fi < \/Mq, we deduce the existence of a 
constant Co > such that, 

For (x, y) G A N , one has s(x) + s(y) + $ > 2S + c (Nk) 2/3 . 

Therefore, by the Calderon-Vaillancourt theorem (and also using Proposi- 



tion 



5.2 in order to regularize the function s(x)), we obtain, 

||e- s /' l ^ J v,i? e ]e- s/ ' l (l-^Ar)(^n)- 1 || 

M\{l-^ N )e- s l h [^ N ,R e ]e- s l h {hD n )- 1 \\ = 0(h c °e- 2S l h ). 

Then, writing, 

(ty N ,R e }u ,i> N u) = {e- s l h {^ N ,Re\e- s l h {e s ' h u e )^ N e s ' h VL), 



and using Proposition 7.1 and Remark 7.2, the result follows. 



□ 



Inserting (8.3) into (8.2), and approaching ipN~u by ipN^-CN, we obtain, 
(8.4) Lmpi = Im (2h 2 (V^ N )Vu CN + h 2 (Aip N )u,ip N u CN ) 



+h{ijj N [ip N , R]i/j N ucN,tpN^CN) + 0(h 



L-N \ e -2S/h_ 



Finally, using Proposition |6.2| (in particular the expression (4.55) of ucn 
in Wn{i) f~l suppV'Tv), we can perform a stationary-phase expansion in 

\SA\ (as in |FLM| . Section 7), and, for L large enough, we obtain, 

L 

Im pi = -h^-TV 2 E hi,mh j+i \ \nh\ m e- 2S ' h + 0(h L / 2 )e- 2S /\ 

j=n 0<m<£<L 

with /n , o,o > 0. In particular, the result for p\ follows. 

The result for pj, j > 2 can be done along the same lines, by using a 



representation of Im pj analogous to (8.1), and by approaching u by a lin- 
ear combination of WKB expressions similar to uq n (where the number of 
terms depends on the asymptotic multiplicity of the resonance: See [HeSj2|, 
Section 10). 
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